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ABSTRACT 

We present the analysis of the light curves of 9 high-magnification single-lens gravitational microlensing 
events with lenses passing over source stars, including OGLE-2004-BLG-254, MOA-2007-BLG-176, MOA- 
2007-BLG-233/OGLE-2007-BLG-302, MOA-2009-BLG-174, MOA-2010-BLG-436, MOA-201 l-BLG-093, 
MOA-2011-BLG-274, OGLE-2011-BLG-0990/MOA-2011-BLG-300, and OGLE-2011-BLG-llOl/MOA- 
201 l-BLG-325. For all events, we measure the linear limb-darkening coefficients of the surface brightness 
profile of source stars by measuring the deviation of the Ught curves near the peak affected by the finite-source 
effect. For 7 events, we measure the Einstein radii and the lens-source relative proper motions. Among them, 5 
events are found to have Einstein radii less than 0.2 mas, making the lenses candidates of very low-mass stars 
or brown dwarfs. For MOA-2011-BLG-274, especially, the small Einstein radius of 0e ^ 0.08 mas combined 
with the short time scale of Ie ^ 2.7 days suggests the possibility that the lens is a free-floating planet. For 
MOA-2009-BLG-174, we measure the lens parallax and thus uniquely determine the physical parameters of 
the lens. We also find that the measured lens mass of 0.84 Mq is consistent with that of a star blended with 
the source, suggesting that the blend is likely to be the lens. Although we find planetary signals for none of 
events, we provide exclusion diagrams showing the confidence levels excluding the existence of a planet as a 
function of the separation and mass ratio. 
Subject headings: gravitational lensing: micro - Galaxy: bulge 



1. INTRODUCTION to the projected lens- source separation by 

When an astronomical object (lens) is closely aligned with 
a background star (source), the light from the source is de- 
flected by the gravity of the lens, resulting in brightening of 

the source star. The magnification of the source flux is related + 2 
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Fig. 1.— Light curve of OGLE-2004-BLG-254. The lower two panels 
show the overall shape of the light curve and residual from the best-fit model. 
The upper two panels show the enlargement of the peak region enclosed by a 
small box in the lower panel. We note that a model light curve varies depend- 
ing on an observed passband due to the chromaticity caused by the finite- 
source effect. The presented model curve is based on the passband of the first 
observatory in the list. However, the residuals of the individual data sets are 
based on the model curves of the coiTesponding passbands. Colors of data 
points are chosen to match those of the labels of observatories where data 
were taken. The two dotted vertical lines in the upper panel represent the 
limb-crossing start/end times. The peak source magnification Ap is given in 
the upper panel. 

where the separation u is expressed in units of the angular Ein- 
stein radius 6^. The Einstein radius is related to the physical 
parameters of the lens system by 



0E = (KM7r,-el) 



1/2. 



TTi-el = AU 



1 



(2) 



where k = 4G/(c^AU) = 8.14 mas Mq~\ M is the mass of the 
lens, TTrei is the relative source-lens parallax, and Dl and Ds 
are the distances to the lens and source star, respectively. The 
relative motion between the source, lens, and observer leads to 
light variation of the source star (lensing ev ent). T he first mi- 
crolensing eve nts w ere detected by lAlcock et al.l (T993) and 
lUdalski et aTl (I1993|) from the exper iments conducted based 
on the proposal of iPaczvnskil (119861) . With the development 
of observational strategy combined with upgraded instrument, 
the detection rate of lensing events has been dramatically in- 
creased from several dozen events per year during the early 
phase of lensing experiments to more than a thousand events 
per year in current experiments. 

The magnification of source star flux increases as the lens 
approaches closer to the source star. For a small fraction of 
events, the lens-source separation is even smaller than the 
source radius and the lens passes over the surface of the source 
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Light curve of MOA-2007-BLG-176. Notations are same as in 



Star These events are of scientific importance due to various 
reasons. 

First, a high-magnification event with a lens passing 
over a source star provides a rare chance to measure the 
brightness profile of a remote star For such an event, 
in which the lens-source separation is comparable to the 
source size near the peak of the event, different parts of 
the source star are magnified by different amounts. The 
resulting lensing light curve deviates from the standard 
form of a point -source event (Witt & Mao 1994; Gould 
ll994tlNemiroff & Wickramasinghe 1994; Alcock et al..l99i) 
and the analysis of the deviation enables to measure 
the limb-darkenin g pro fi le of the lens e d star ( Witli 119951: 
Loeb & Sasselo\l Il99l IValls-GabaudI [T99l [Brvce et al.1 
2002: Heyrovsk?! I2003h . With the same principle, it is also 
possible to study irregular sur face structures such as spots 
(Hevrovskv & Sasselov 2000: Han et al.ll2000l: iHendry et all 
12002; .Rattenburv et al.. .2002 ) . 

Second, it is possible to measure the Einstein radius of the 
lens and the relative lens-source proper motion. The light 
curve at the moment of the entrance (exit) of the lens into 
(from) the source surface exhibits inflection of the curvature. 
The duration of the passage over the source as measured by 
the interval between the entrance and exit of the lens over the 
surface of the source star is 



AfT = 2\/ p*--uo^ t-E, 



(3) 



where p^, is the source radius in units of 9^ (normalized source 
radius), mq is the lens-source separation normalized by 9^, at 
the moment of the closest approach (impact parameter), and 
is the time scale for the lens to transit 9e (Einstein time 
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SOURCE-TRANSIT LENSING EVENTS 



TABLE 1 

Events With Lenses Passing Over Source Stars 



event 



reference 



OGLE-2()04-BLG-254 

OGLE-2()04-BLG-482 

MOA-2006-BLG- 1 30/OGLE-2006-BLG-437 

OGLE-2()07-BLG-050/MOA-2007-BLG- 1 03 

MOA-2007-BLG-176 

OGLE-2()07-BLG-224/MOA-2007-BLG- 1 63 

MOA-20()7-BLG-233/OGLE-2007-BLG-302 

OGLE-2()08-BLG-279/MOA-2008-BLG-225 

OGLE-2008-BLG-290/MOA-2008-BLG-241 

MOA-20()9-BLG-174 

MOA-2009-BLG-41 1 

MOA-2010-BLG-311 

MOA-2010-BLG-436 

MOA-2010-BLG-523 

MOA-2011-BLG-093 

MOA-2011-BLG-274 

OGLE-2()l 1-BLG-()990/MOA-201 l-BLG-300 
OGLE-2() 1 1 -BLG- 1 1 l/MOA-20 1 1 -BLG-325 



Cassan et al. (2006'! / this work 
Zub etal. (2011) 

Baudrv et al. (201 1) / under analysis 
Batista et al. (2009) 
this work 

Gould et al. (2009) 
this work 
Yee et al. (2009) 
Fouque et al. (2010) 
this work 

J'ouaue et al. (20 11) / under analysis 
^lung et al. (201 1) / under analysis 
this work 

^ouldetak QOl^ / under analysis 
this work 
this work 
this work 
this work 



scale). The impact parameter and the Einstein time scale are 
measured from the overall shape of the light curve and the 
duration of the event. With the known uq and Ie combined 
with the measured duration of passage over the source, the 
normalized source radius is measured from the relation (O. 
With the additional information about the angular source size, 
6**, then the Einstein radius and the lens-source proper motion 
are measured as Be = ^^,/p^, and /i = Oe/Ie, respectively. For 
general lensing events, the Einstein time scale is the only mea- 
surable quantity related to the physical parameters of the lens. 
However, the time scale results from the combination of 3 
physical parameters of the mass of the lens, M, the distance to 
the lens, De, and the lens-source transverse speed, v, and thus 
the information about the lens is highly degenerate. The Ein- 
stein radius, on the other hand, does not depend on v and thus 
the physical parameters of the lens can be better constrained. 
For a fraction of events with long time scales, it is possible 
to additionally measure the lens parallax, tte = t^kx/Oe, from 
the deviation of the light curve induced by the orbital motion 
of the Earth around the Sun. With the Einstein radius and the 
lens parallax measured, the physical parameters of the lens 
are uniquely determined (Gould 1997|). 

Third, high-magnification events are sensitive to planetary 
companions of lenses. This is because a planet induces a small 
caustic near the primary lens and a high-magnification event 
resulting from the source trajectory passing close to the pri- 
mary has a high chance to produce signals indicating the ex- 
istence of the planet dGriest & Safizadeh 1998). For an event 
with a lens passing over a source star , the planeta ry signal is 
weakened by the finite-source effect (iBennett & R hie 1996j). 
Nevertheless, two of the microlensing planets were discov- 
ered through this channel: MOA- 2007-BLG-4 0Q (Dong et all 
[2009') and MOA-2008-BLG-310 (J anczak et al. 2010). 

Fourth, high-magnification events provide a chance to spec- 
troscopically study remote Galactic bulge stars. Most stars 
in the Galactic bulge are too faint for spectroscopic observa- 
tions even with large telescopes. However, enhanced bright- 
ness of lensed stars of high-magnification events allows spec- 
troscopic observation possible, enabling population sti idy of 
Galactic bulge star s (Johnson et al. 2008; Be nsbv et^l2009l 
1201 UlCohen et al.ll2C)09; Epstein et al. .2010). 

In this work, we present integrated results of analysis for 
14 high-magnification events with lenses passing over source 
stars that have been detected since 2004. Among them, 8 



events were newly analyzed and one event was reanalyzed 
with additional data. 

2. EVENT SELECTION 

The sample of events in our analysis is selected under the 
definition of a single-lens event where the lens-source sepa- 
ration at the time of the peak magnification is less than the 
radius of the source star, i.e. uq < and thus the lens passes 
over the surface of the source star To obtain a sample of 
events, we begin with searching for high-magnification events 
that have been detected since 2004. Events with lenses pass- 
ing over source stars can be usually distinguished by the char- 
acteristic features of their light curves near the peak. These 
features are the inflection of the curvature at the moment when 
the finite source first touches and completely leaves the lens 
and the round shape of the light curve during the passage of 
the lens over the source. To be more objective than visual 
inspection, we conduct modeling of all high-magnification 
events with peak magnifications Ap > 10 to judge the qual- 
ification of events. From these searches, we find that there 
exist 18 such events. Among them, analysis results of 12 
events were not published before. We learn that 4 unpub- 
hshed events MOA -2006-BLG- 1 30/OGLE-2006-B LG-437 
(iBaud rv et all I20TTI) . MOA- 2009-BLG-411 (.Fouqueetall 
201 1), MOA-201 0-BLG-523 dGould et al.|[201 li and MOA- 

2010- BLG-311 (lHungetal.1 120111) are under analysis by 
other researchers and thus exclude them in our analy- 
sis. We note that there exist 4 known source-crossing 
events detected befo re 2004, including M ACHO Aler t 95-30 
(lAlcocketal.lll997h . O GLE sc26 2218 dSmitheLall [2003), 
OGLE-2003-BLG-238 ("Jiang et al.ll2004 l). and OGLE-2003- 
BLG-262 (Yooetal. 2004)' We also note that MOA- 
2007-B LG-400 (Dong et al. 2009) and MOA-2008-BLG-310 
(I Janczak et al. 2010) exhibit characteristic features of source- 
crossing single-lens events but we exclude them in the sample 
because the lenses of the events turned out to have planetary 
companions. 

In this work, we conduct analyses of 9 events. Among 
them, 8 events are newly analyzed in this work. These events 
include MOA-2007-BLG-176, MOA-2007-BLG-233/OGLE- 
2007-BLG-302, MOA-2009-BLG-174, MOA-2010-BLG- 
436, MOA-2011-BLG-093, MOA-2011-BLG-274, OGLE- 

201 1- BLG-0990/MOA-201 l-BLG-300, and OGLE-2011- 
BLG-llOl/MOA-BLG-201 1-325. For OGLE-2004-BLG- 
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TABLE 2 
Observatories 



event 
(RA,DEC)j2000 


MOA 


OGLE 


^tPUN 


PLANET 


RoboNet 
/MiNDSTEp 


OGLE-2004-BLG-254 
(17''56"'36".20,-32°33'0I".80) 




LCOi 
(377) 


CTIOi.v (39/5) 
FCOn (129) 


Boydeni (74) 
Canopusi (59) 
SAAOi (112) 




MOA-2()07-BLG-I76 
(18''05"'00MI ,-25°47'03".69) 


Mt.JohnR 
(1388) 




AucklandR (68) 
CTIOi V (41/4) 
FCOn (33) 
Lemmoni (66) 
VLOn (129) 


Canopusi (26) 
Stewardi (4) 


Danishi (2) 


MOA-2007-BLG-233 

/OGLE-2007-BLG-302 

(17''54"'I4".86-31°I1'02".65) 


Mt.JohnR 
(645) 


LCOi 
(628) 


CTIOi V (80/5) 
FCOn (23) 
Lemmoni (19) 
SSOn (80) 


Canopusi V (60/5) 
Perthi (23) 
SAAOi (80) 


Danishi (125) 


MOA-2009-BLG-I74 

( 1 8''02"07''.60 -3 1 ° 25'24".20) 


Mt.JohnR 
(2189) 




BronbergN (147) 
CAOn (111) 
CraigieN (130) 
CTIOi.v (286/7) 
KumeuN (90) 
PossumN (60) 


Canopusi (40) 


LTr (7) 


MOA-2()I0-BLG-436 

(18^03"21!68,-27°38'10".74) 


Mt.JohnR 
(2581) 






SAAOi.v (14/3) 


FTSi (3) 


MOA-2()ll-BLG-()93 

(17''46"'I7".83,-34°20'24".76) 


Mt.JohnR 
(2247) 


LCOi 
(292) 


CTlOi.v (76/21) 
PESTn (124) 


Canopusi (254) 


FTNi (3) 
FTSi (19) 


MOA-2()II-BLG-274 

(17h54"i42";34-28°54'59".26) 


Mt.JohnR 

(3447) 


LCOi 

(76) 


AucklandR (53) 
CTIOi (4) 
FCOn (16) 
KumeuR (49) 
PESTn (15) 






OGLE-20II-BLG-0990 
/MOA-20I1-BLG-300 
(I7''5I"30";29 -30° I7'47".60) 


Mt.JohnR 
(1708) 


LCOi 
(3434) 


OPDi (275) 
PossumR (23) 


Canopusi (10) 
SAAOi.v (95/6) 




OGLE-20II-BLG-II0I 

/MOA-20I1-BLG-325 

(18''03""31';62,-26°20'39".50) 


Mt.JohnR 
(609) 


LCOi 
(192) 


AucklandR (60) 
CTIOi.v (126/12) 
PossumR (24) 
SSOn (107) 
VLOn (113) 


Canopusi (98) 


FTNi (65) 
FTSi (145) 
LTi (27) 



Note. — Mt. John: Mt. John Observatory, New Zealand; LCO: Las Campanas Observatory, Chile; 
Auckland: Auckland Observatory, New Zealand; Bronberg: Bronberg Observatory, South Africa; 
CAO: CAO San Pedro Observatory, Chile; Cragie: Craigie Observatory, Australia; CTIO: Cerro Tololo 
Inter-American Observatory, Chile; FCO: Farm Cove Observatory, New Zealand; Kumeu: Kumeu Ob- 
servatory, New Zealand; Lemmon: Mt Lemmon Observatory, Arizona, USA; OPD: Observatorio do 
Pico dos Dias, Brazil; PEST: Perth Exoplanet Survey Telescope, Australia; Possum: Possum Observa- 
tory, New Zealand; SSO: Southern Stars Observatory, French Polynesia; VLO: Vintage Lane Observa- 
tory, New Zealand; Boyden: Boyden Observatory, South Africa; Canopus: Canopus Hill Observatory, 
Tasmania, Australia; Perth: Perth Observatory, Australia; SAAO: South African Astronomical Ob- 
servatory, South Africa; Steward: Steward Observatory, Arizona, USA; FTN: Faulkes North, Hawaii; 
FTS: Faulkes South, Australia; LT: Liverpool Telescope, La Palma, Spain; Danish: Danish Telescope, 
European Southern Observatory, La Silla, Chile. The subscription after each observatory represents 
the filter used for observation and the value in parenthesis is the number of data points. The filter "N" 
denotes that no filter is used. 



254, which was analyzed before by iCassan et alj (12006 ^. we 
conduct additional analysis by adding more data sets taken 
from CTIO and FCO. Q In Table [T] we summarize the status 
of analysis for the total 1 8 events that have been detected since 
2004. 

3. OBSERVATION 

For almost all events analyzed in this work, the source- 
crossing part of the light curve was densely covered. This was 
possible due to the coordinated work of survey and follow- 

' Besides the data sets Hsted in Table [2] there exists an additional data 
set taken by using the Danish telescope. However, we do not use these data 
because it has been shown by Heyrovsky 1 2008) that the large scatter of the 
data results in poor measurement of tensing parameters including the limb- 
darkening coefficient. 



up observations. Survey groups issued alerts of events. For 
a fraction of the events with high-magnifications, additional 
alerts were issued. In other cases, follow-up teams issued 
high-magnification alerts independently. The peak time of a 
high-magnification event was predicted by real-time model- 
ing based on the rising part of the light curve. Finally, the 
peak was densely covered by many telescopes that were pre- 
pared for follow-up observations at the predicted time of the 
peak. For MOA-2010-BLG-436, the rising part of the light 
curve was not covered by survey observations due to the short 
time scale of the event and thus no alert was issued. Nev- 
ertheless, the event was positioned in a high frequency field 
of the MOA survey and thus the peak was covered densely 
enough to be confirmed as an event with the lens passing over 
the source. 
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M0A-Z007-BLG-a33/OGLE-2007-BLG-30Z 




4370 4280 4290 4300 4310 

HJD-2450000 



Fig. 3.— Light curve of MOA-2007-BLG-233/OGLE-2007-BLG-302. 
Notations are same as in Fig. 1 

Table |2] shows the observatories of the telescopes that were 
used for observations of the individual events along with the 
observed passbands (marked as subscripts after the observa- 
tory names) and the numbers of data points (values in paren- 
theses). Also marked are the coordinates (RA,DEC) of the 
events. Survey observations were conducted by MOA and 
OGLE groups using the 1.8 m telescope at Mt. John Ob- 
servatory in New Zealand and the 1.3m Warsaw University 
telescope at Las Campanas Observatory in Chile, respectively. 
Follow-up observations were carried out by /iFUN, PLANET, 
RoboNet, and MiNDSTEp groups using 22 telescopes lo- 
cated in 8 different countries. These telescopes include 1.3 
m SMARTS CTIO, 0.4 m CAO in Chile, 0.4 m Auckland, 0.4 
m FCO, 0.4 m Possum, 0.4 m Kumeu, 0.4 m VLO in New 
Zealand, 1.0 m Lemmon in Arizona, USA, 0.4 m Bronberg in 
South Africa, 0.6 m Pico dos Dias in Brazil, 0.25 m Craigie, 
0.3 m PEST in AustraUa, 0.28m SSO in French Polynesia, 1.0 
m SAAO, 1.5 m Boyden in South Africa, 1.0 m Canopus, 0.6 
m Perth in Australia, 1.5 m Steward in Arizona, USA, 2.0 m 
FTN in Hawaii, USA, 2.0 m FTS in AustraHa, 2.0 m LT in La 
Palma, Spain, and 1.54 m Danish in La Silla, Chile. 

Reduction of data was conducted by using photometry 
codes that were developed by the individual groups. The 
MOA and O GLE data were re duce d by photomet ry codes 
developed by lBondetall (2001) and lUdalski l (2003), respec- 
tively, which are based on Difference Image Analysis method 
dAlard & Lupton 1991. The mFUN data were processed us- 
ing a DoPHOT pipeline (I Schechteii[T993D. For PLANET a nd 
MiNDSTEp data, a pySIS pipeline dAlbrow et alj 119981) is 
used. For RoboNet data, a DanDIA pipeUne ( lBramichll2008l) 
is used. 



M0A-Z009-BLG-174 




4940 4950 4960 4970 4980 



HJD-2450000 

Fig. 4. — Light curve of MOA-2009-BLG-174. Notations are same as in 
Fig.l 

The error bars estimated from different observatories are 
rescaled so that pei" degree of freedom becomes unity for 
the data set of each observatory where is computed based 
on the best-fit model. According to this simple scheme, how- 
ever, we find a systematic tendency for some data sets that 
error bars near the peak of a light curve are overestimated. 
We find that this is caused by the inclusion of redundant data 
at the baseline in error normalization. In this case, the data at 
the baseline greatly outnumber accurate data points near the 
peak and thus error-bar normalization is mostly dominated by 
the baseline data. To minimize this systematics, we restrict 
the range of data for error normalization not to be too wide 
so that error estimation is not dominated by data at the base- 
line, but not to be too narrow so that lensing parameters can 
be measured accurately. For the final data set used for model- 
ing, we eliminate data points lying beyond 3(7 from the best-fit 
model. 

In Figure [T]|9] we present the light curves of the individual 
events. In each figure, the lower two panels show the over- 
all shape of the light curve and residual and the upper two 
panels show the enlargement of the peak region of the light 
curve and residual. For each figure, we mark the moments of 
the lens'entrance and exit of the source by two dotted vertical 
lines. Also marked is the peak source magnification. We note 
that the same color of data points is used for each observa- 
tory throughout the light curves and colors of data points are 
chosen to match those of labels of observatories. We note that 
the magnitude scale corresponds to one of the observatories in 
the list, while data from the other observatories have adjusted 
blends and are vertically shifted to match the first light curve. 
The choice of reference is based on data from survey obser- 
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Fig. 5. — Light curve of MOA-201()-BLG-436. Notations are same as in 
Fig.l 

vation, i.e. OGLE and MOA data. If both OGLE and MOA 
data are available, the OGLE data is used for reference. 

4. MODELING 

Modeling the light curve of each event is conducted by 
searching for a set of lensing parameters that best describes 
the observed light curve. For all events, the light curves ap- 
pear to have a standard form except the peak region and thus 
we start with a simple single-lens modeling. The light curve 
of a standard single-lensing event is characterized by 3 param- 
eters, that are the time of the closest lens-source approach, fo, 
the lens-source separation at that moment, mq, and the Ein- 
stein time scale, Ie- Based on the initial solution, we refine 
the solution by considering additional second-order effects. 

To precisely describe the peak region of the light curve of 
an event with a lens passing over a source, additional parame- 
ters are needed to describe the deviation caused by the finite- 
source effect. To the first order approximation, the finite- 
source effect is described by the normalized source radius, p*. 
For more refined description of the deviation, additional pa- 
rameters of the limb-darkening coefficients, ux, are needed to 
account for the variation of the deviation caused by the bright- 
ness profile of the source star surface. With the coefficients, 
the limb-darkening profile is modeled by the standard linear 
law 

I = Io[l-ux(l-cos(b)], (4) 

where /q is the intensity of light at the center of the stellar disk 
and (j) is the angle between the normal to the stellar surface 
and the line of sight toward the observer. 

For an event with a time scale comparable to the orbital pe- 
riod of the Earth, the position of the observer changes by the 
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5660 5670 5680 5690 5700 

HJD-e450000 

Fig. 6. — Light curve of MOA-201 l-BLG-093. Notations are same as in 
Fig.l 

Earth's orbital motion during the event and the resulting light 
curve deviates from a symmetric standard form. This parallax 
effect is described by 2 parameters of tte,n and tte.b, that rep- 
resent the two components of the lens parallax vector tte pro- 
jected on the sky in the north and east equatorial coordinates, 
respectively. The direction of the parallax vector corresponds 
to the lens-source relative motion in the frame of the Earth. 
The size of the parallax vector corresponds to the ratio of the 
Earth's orbit, i.e. 1 AU, to the Einstein radius projected on the 
observer's plane. 

For a high-magnification event, the peak of the light curve 
can exhibit additional deviations if the lens has a companion. 
For a planetary companion located near the Einstein ring of 
the primary lens or a binary companion with a separation from 
the primary substantially smaller or larger than the Einstein 
radius, a small caustic is induced near the location of the pri- 
mary lens. Then, the source trajectory of a high-magnification 
event passing close to the primary lens has a high chance to 
approach the caustic, resulting in a perturbation near the peak 
of the light curve. Description of the perturbation induced 
by a lens companion requires 3 additional parameters of the 
mass ratio between the binary lens components, q, the pro- 
jected separation in units of the Einstein radius, s, and the 
angle between the source trajectory and the binary axis, a. 

For each event, we search for a solution of the best-fit 
lensing parameters by minimizing in the parameter space. 
For the minimization, we use a Markov Chain Monte 
Carlo (MCMC) method. We comp ute finite magnifications by 
using the ray -shoo ting technique jSchneider & Weiss! 119861 : 
iKavser et al.T[T986t IWambsganssl ll997h ~ In this method, rays 
are uniformly shot from the image plane, bent according to 
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Light curve of MOA-201 l-BLG-274. Notations are same as in 



the lens equation, and land on the source plane. Then, a finite 
magnification is computed by comparing the number densi- 
ties of rays on the image and source planes. Precise com- 
putation of finite magnifications by using this numerical tech- 
nique requires a large number of rays and thus demands heavy 
computation. To minimize computation, we limit finite- 
magnification computation by using the ray-shooting method 
only when the lens is close to the source. Once a solution of 
the parameters is found, we estimate the uncertainties of the 
individual parameters based on the chain of solutions obtained 
from MCMC runs. 

5. RESULT 

In Table [3] we present the lensing parameters of the best- 
fit solutions of the individual events determined from model- 
ing. To provide integrated results of events with lenses pass- 
ing over source stars, we also provide solutions of events that 
were previously analyzed. For OGLE-2004-BLG-254, we 
provide both solutions of the previous analysis and this work 
for comparison. 

For all events analyzed in this work, we are able to mea- 
sure the limb-darkening coefficients of source stars. In Ta- 
ble |4] we present the measured limb-darkening coefficients. 
We measure the coefficients corresponding to the individual 
data sets covering the peak of each light curve instead of the 
individual passbands. This is because the characteristics of 
filters used for different telescopes are different from one an- 
other even though they are denoted by a single representa- 
tive band and thus joint fitting of data measured in different 
filter systems may result in erroneous measurement of limb- 
darkening coefficients tFouque et al...2010.) . To compare with 




5760 
HJD-2450000 



Fig. 8.— Liglit curve of OGLE-201 1-BLG-0990/MOA-2011-BLG-300. 
Notations are same as in Fig. 1 

theoretical values, we also provide values of coefficients pre- 
dicted by Claret (2000) for the Bessell V, R, and / filters. 
Also provided are the source types and the adopted values 
of log^ and T^ff where the typical uncertainties of the sur- 
face gravity and the effective temperature are A(logg) = 0.5 
and ATeff = 250 K, respectively. We adopt a solar metallicity. 
We note that the measured coefficients are generally in good 
agreement with theoretical values, Mth- From the table, it is 
found that for 23 out of the total 29 measurements the mea- 
sured coefficients are within 20% range of the fractional dif- 
ference as measured by /„ = {u-u±)/utif The cases of large 
differences with /„ > 20% include M/(Canopus) for OGLE- 
2004-BLG-254, mr(MOA), MR(Auckland), and M/(Lemmon) 
for MOA-2007-BLG-176, m,(CTIO) for MOA-2009-BLG- 
174, and m,(CTIO) for OGLE-201 1-BLG-l lOl/MOA-201 1- 
BLG-325. From the inspection of the individual data points 
on the light curves, we find the major reasons for the large dif- 
ferences between the measured and theoretical values are due 
to poor coverage [M/(Canopus) for OGLE-2004-BLG-254, 
Mfi(MOA, Auckland) and M/(Lemmon) for MOA-2007-BLG- 
176, M/(CTIO) for MOA-2009-BLG-174] or poor data qual- 
ity [m/(CTIO) for OGLE-2011-BLG-1101/MOA-2011-BLG- 
325] . Other possible reason for differences from the predicted 
values include differences of individual filters from the stan- 
dard Bessel filters, as well as differences in the method to 
compute the theoretical values (Heyrovsky 2003). 

The source type of each event is determined based on the 
location of the source in the color-magnitude diagram (CMD) 
of stars in the same field. CMDs are obtained from CTIO im- 
ages taken in V and / bands. To locate the lensed star in the 
CMD, it is required to measure the fraction of blended fight in 
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TABLE 3 
Best-fit Parameters 



event 


xVdof 


(0 (HJD-2450000) 


"0 


tE (days) 


P* 




OGLE-2004-BLG-254 


1326/593 


3166.8194 


0.0046 


13.23 


0.0400 


— 






±0.0002 


±0.0008 


±0.05 


±0.0002 


- 


OGLE-2004-BLG-482 


756.9/693 


3235.7816 


0.000 


9.61 


0.1309 


- 






±0.0007 


±0.002 


±0.02 


±0.0005 


- 


OGLE-2007-BLG-050 


1760.5/1745 


4221.9726 


0.002 


68.09 


0.0045 


0.12 


/MOA-2007-BLG- 1 03 




±0.0001 


±0.000 


±0.66 


±0.0001 


±0.03 


OGLE-2007-BLG-224 


- 


- 


0.00029 


6.91 


0.0009 


1.97 


/MOA-2007-BLG-163 




- 


- 


±0.13 


±0.0002 


±0.13 


OGLE-2008-BLG-279 


- 


4617.34787 


0.00066 


106.0 


0.00068 


0.15 


/MOA-2008-BLG-225 




±0.00008 


±0.00005 


±0.9 


±0.00006 


±0.02 


OGLE-2008-BLG-290 


2317.7/2015 


4632.56037 


0.00276 


16.36 


0.0220 


- 


/MOA-2008-BLG-241 




±0.00027 


±0.0002 


±0.08 


±0.0001 


- 


OGLE-2004-BLG-254 


785.2/784 


3166.823 


0.0111 


12.84 


0.0418 


- 






±0.001 


±0.0004 


±0.09 


±0.0004 


- 


MOA-2007-BLG-176 


1756.0/1747 


4245.056 


0.0363 


8.13 


0.0590 


- 






±0.001 


±0.0005 


±0.07 


±0.0006 


- 


MOA-2()07-BLG-233 


1779.4/1757 


4289.269 


0.0060 


15.90 


0.0364 


- 


/OGLE-20()7-BLG-302 




±0.001 


±0.0002 


±0.05 


±0.0001 


- 


MOA-2009-BLG-174 


2816.5/3051 


4963.816 


0.0005 


64.99 


0.0020 


0.06 






±0.001 


±0.0001 


±0.61 


±0.0001 


+0.07-0.02 


MOA-2010-BLG-436 


2599.4/2593 


5395.791 


0.0002 


12.78 


0.0041 








±0.001 


±0.0002 


±1.08 


±0.0003 




MOA-2011-BLG-093 


3038.0/3024 


5678.555 


0.0292 


14.97 


0.0538 








±0.001 


±0.0002 


±0.05 


±0.0002 




MOA-2011-BLG-274 


3657.7/3649 


5742.005 


0.0029 


2.65 


0.0129 








±0.001 


±0.0001 


±0.06 


±0.0003 




OGLE-2011-BLG-0990 


5551.6/5540 


5758.691 


0.0151 


6.70 


0.0199 




/MOA-2011-BLG-300 




±0.001 


±0.0004 


±0.07 


±0.0003 




OGLE-2011-BLG-llOl 


1562.6/1562 


5823.574 


0.0485 


29.06 


0.0979 




/MOA-2011-BLG-325 




±0.002 


±0.0005 


±0.11 


±0.0006 





Note. — The parameters of the first 6 events are adopted from previously analyses and those of 
the other 9 events are determined in this work. For OGLE-2004-BLG-254, the event was reanalyzed by 
adding more data sets. The references of the previous analyses are presented in Table[T] 



the observed light curve. This is done by including a blending 
parameter in the process of light curve modeling. For MOA- 
2011-BLG-274, a CMD taken from CTIO is available but im- 
ages were taken after the event and thus we could not deter- 
mine the source color and magnitude by the usual method. 
Instead we employ the method of Gould et alj (1201 01). In this 
method, we first measure the source instrumental magnitudes 
by fitting the OGLE (/qgle) and PEST (unfiltered, A^pest) 
data to the light curve model. We then align each of these 
data sets to CTIO {V /I) using comparison stars, which effec- 
tively transforms A^pest/^ogle to (y//)cTio- In figure [TO] we 
present the CMDs of stars in the fields of the individual events 
and the locations of source starsQ For MOA-2010-BLG-436 
and OGLE-201 1-BLG-0990/MOA-201 l-BLG-300, there ex- 
ists SAAO data taken in / and V bands, but the number and 
quality of V-band data are not numerous and good enough to 
specify the source type. 

In Table |5] we present the measured Einstein radii. The 
Einstein radius of each event is determined from the angu- 
lar source radius, 0*, and the normalized source radius, p^, 
as 9e = 9*/ Pi,- The normalized source radius is measured 
from modeling. To measure the angular source radius, we 
use the method of Yoo et al.. (,2004) , where the de-reddened 
V-I color is measured from the location of the source in 

^ We note that high-resolution spectra are available for some events 
with lenses passing over source stars. These events are OGLE-2004- 
BLG-254 (Cassanetal. 2006), OGLE-2004 -BLG-482 <Zub et al. 2011), 
OGLE-2007-BLG-050/MOA-2007-BLG-103 gohnson et aiTBoilD. MOA- 
2009-BLG- 1 74, MOA-2Q10-BLG -31 1, MOA-2010-BLG -523 ffiensbv et alj 
IMIll) . and MOA-201 l-BLG-093 JMcGregor et aU201 in . For those who are 
more interested in the source stars of these events, see the related references. 



th e CMD, V-I is conv erted into V -K using the relation 
of iBessel & Brett! ( 119881) , and then the angular source radius 
is inferred from the V-K col or and the surface brightness 
relation given by Ker vella et a l. (2004). In this process, we 
use the centroid of bulge clump giants as a reference for the 
calibration of the color and brightness of a source under the 
assumption that the source and clump giants experience the 
same amount of extinction and reddening. We note that no 
CMD is available for MOA-201 0-BLG-436 and OGLE-201 1- 
BLG-0990/MOA-2011-BLG-300 and thus the Einstein radius 
is not provided. Also provided in Table|5]are the relative lens- 
source proper motions as measured by p = dE/tE- 

We note that the measured Einstein radii of some events 
are substantially smaller than a typical value. These events 
include OGLE-2004-BLG-254 (Oe - 0.14 mas), MOA-2007- 
BLG-176 (- 0.14 mas), MOA-2007-BLG-233/OGLE-2007- 
BLG-302 (- 0.17 mas), MOA-2011-BLG-093 (- 0.07 mas), 
and MOA-2011-BLG-274 0.08 mas). The lens mass and 
distance are related to the Einstein radius by 

M = 0.019mJ^V^V^^V (5) 
°V8kpcy V^s-£)l/ V0-14masy 

Hence, the small 9e of these events implies that lenses are ei- 
ther very close to the source or very low-mass objects. Most 
of these events have proper motions that are typical of bulge 
lenses (2-7 masyr"') and so may be quite close to the source 
(see Table 5). But MOA-2011-BLG-274 has a substantially 
higher proper motion, /i ^ 11 mas yr"' . It is therefore a 
good candidate for a sub-stellar object or even a free-floating 
planet (,Sumi et al.,,201 1,) . Because of its high proper motion. 
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TABLE 4 
Source Parameters 



event 


source type 


limb-darkening coefficients 




(log Tiff) uv 


I'R 


«/ 


OGLE-2004-BLG-254 


Kill 


0.70±0.05 


0.55±0.05 


OGLE-2004-BLG-482 


Mill 


0.88±0.02 


0.71±0.01 


OGLE-2007-BLG-050 


subgiant - 






/MOA-2007-BLG-103 










r V — 


- 


- 


/MUA-zUU /-oLLr- ioJ 








OGLE-2008-BLG-279 


GV 


- 


- 


/iViUrt-ZUUooljO-ZZ J 










JS.111 U. / / itU.Ui 


0.62±0.07 


0.55±0.01 


/MOA-2008-BLG-241 








OGLE-2004-BLG-254 


Kill 

/"o n /i7^n i^t 
(z.U, 4/jU t^) 


- 


0.70±0.07 (OGLE) 
0.56±0.10(CTIO) 
0.69±0.10(Boyden) 
0.78±0.09 (Canopus) 
0.55±0.06 (SAAO) 
0.61 (Claret 2000) 


MOA-2007-BLG-176 


Kill 


0.53±0.04 (MOA) 


0.50±0.05 (CTIO) 




(2.0, 4500 K) 


0.51 ±0.05 (Auckland) 
0.73 (Claret 2000) 


0.44±0.06 (Lemmon) 
0.63 (Claret 2000) 


MOA-2007-BLG-233 


GUI 


0.56±0.02 (MOA) 


0.53±0.04 (OGLE) 


/OGLE-2007-BLG-302 


(2.5, 5000 K) 


0.68 (Claret 2000) 


0.56±0.02 (CTIO) 
0.49±0.02 (Canopus) 
0.59 (Claret 2000) 


MOA-2009-BLG-174 


FV 

(4.5, 6750 K) 




0.33±0.02 (CTIO) 
0.46 (Claret 2000) 


MOA-2010-BLG-436 




0.52±0.10(MOA) 





MOA-2011-BLG-093 



MOA-2011-BLG-274 



OGLE-2011-BLG-0990 
/MOA-2011-BLG-300 
OGLE-2011-BLG-llOl 
/MOA-2011-BLG-325 



GUI 

(3.0, 5500 K) 



GV 

(4.0, 6000 K) 



Kill 

(2.0, 4250 K) 



0.69±0.05 (CTIO) 
0.70 (Claret 2000) 



0.89±0.14(CTIO) 
0.83 (Claret 2000) 



0.55±0.04 (MOA) 
0.63 (Claret 2000) 



0.48±0.02 (Kumeu) 
0.51±0.01 (Auckland) 
0.S9 (Claret 2000) 



0.77±0.08 (MOA) 
0.76 (Claret 2000) 



0.51 ±0.10 (OGLE) 
0.58±0.04 (CTIO) 
0.51 ±0.03 (Canopus) 
0.54 (Claret 2000) 



0.56±0.04 (OPD) 

0.74±0.07 (OGLE) 
0.81 ±0.07 (CTIO) 
0.77±0.06 (Canopus) 
0.78±0.05 (FTS) 
0.65 (Claret 2000) 



Note. — The parameters of the first 6 events are adopted from previously analyses and those of 
the other 9 events are analyzed in this work. For OGLE-2004-BLG-254, the event was reanalyzed by 
adding more data sets. The references of the previous analyses are presented in Table [T] The limb- 
darkening coefficients, u\, are presented for the in dividual data sets used for u\ measurements and they 
are compared with theoretical values predicted bv lClarell ( 120001) . Also presented are the adopted values 
of logg and Tcff. The unit of the stellar surface gravity is cm/s^ 



it should be possible to detect the lens within a few years us- 
ing high-resolution infrared imaging, provided it is luminous. 
In this case a null result would confirm its substellar nature. 

For MOA-2009-BLG-174, the lens parallax is measured 
with Ax^ ^ 16.2. The measured parallax parameters are 

tte II = -0.049 ± 0.006 ; tte.x = 0.038 ± 0.065, (6) 

where tt^ || and tte.± are the components of the lens par- 
allax vector that are parallel with and perpendicular to the 
projected position of the Sun. These values correspond to 
the standard parallax components of {TrE.N,'^E,E) = (0.025 ± 
0.052,-0.057 ±0.028). In Figure fTTI we present contours of 
in the space of the parallax parameters. Combined with the 
measured Einstein radius, the physical parameters of the lens 
are uniquely determined as 



and 



AU 



:6.39±l.llkpc, 



(8) 



M= — ^ =0.84 ±0.37 Mq, 

KTTe 



(7) 



respectively. We find that the measured lens mass is consistent 
with the de-reddened color of blended light (V-I)qi, ~ 1.4, 
which approximately corresponds to the color of an early K- 
type main-sequence star with a mass equivalent to the esti- 
mated lens mass, suggesting that the blend is very likely to be 
the lens. We mark the position of the blend in the correspond- 
ing CMD in Figure [To] 

A high-magnification event is an important target for planet 
search due to its high efficiency to planetary perturbations. 
Unfortunately, we find no statistically significant deviations 
from the single-lens fit for any of the events analyzed in this 
work. However, it is still possible to place limits on the range 
of the planetary separation and mass ratio. For this purpose, 
we construct so-call "exclusion diagrams" which show the 
confidence levels of excluding the existence of a planet as 
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TABLE 5 
Physical Lens Parameters 



event 


Oe (mas) 


fi (mas yr ^ ) 


M(Mq) 


DL(kpc) 


WLii^j:l~ZUU4-i3Hj-2j4 


n 1 1 /I 
U. 1 14 


"2 1 






r~\r^T c onri/i dt /"* /iot 


U.4 


1 A 






UOLiil-ZUU /-dLO-UjU/MUA-zUU /-oLLj- IUj 


U.4fSdiU.Ul 


Z.ujiLU.Uo 


U.jUitU. 14 


J. JltU.4 


OGLE-2007-BLG-224/MOA-2007-BLG- 1 63 


0.91ib0.04 


48±2 


0.056±0.004 


0.53±0.04 


OGLE-2008-BLG-279/MOA-20()8-BLG-225 


0.81 ±0.07 


2.7±0.2 


0.64±0.10 


4.0±0.6 


OGLE-2008-BLG-290/MOA-20()8-BLG-241 


0.30ib0.02 


6.7±0.4 






OGLE-2004-BLG-254 


0.14±0.01 


4.()6±0.35 






MOA-2007-BLG-176 


O.UibO.Ol 


6.21±0.54 






MOA-2007-BLG-233/OGLE-2007-BLG-302 


0.17ib0.01 


3.81±0.33 






MOA-2009-BLG-174 


0.43ib0.04 


2.40±0.24 


0.84±0.37 


6.39±L11 


MOA-2010-BLG-436 










MOA-2011-BLG-093 


O.OTibO.Ol 


L80±0.16 






MOA-2011-BLG-274 


O.OSibO.Ol 


1L18±0.97 






OGLE-201 1-BLG-0990/MOA-201 l-BLG-300 










OGLE-201 1-BLG-l lOl/MOA-201 l-BLG-325 


0.24ib0.02 


2.99±0.26 







Note. — The parameters of the first 6 events are adopted from previously analyses and 
those of the other 9 events are analyzed in this work. For OGLE-2004-BLG-254, the event was 
reanalyzed by adding more data sets. The references of the previous analyses are presented in 
Tablelll 




5820 5B40 
HJD-245000D 



Fig. 9.— Light curve of OGLE-2011-BLG-1101/MOA-2011-BLG-325. 
Notations ai'e same as in Fig. 1 



a function of the normalized star-planet separation and the 
planet/star mass ratio. We construct diagrams by adopting 
iGaudi & Sacketl! (2000) method. In this method, binary mod- 
els are fitted to observed data with the 3 binary parameters 
(s, q, a) are held fixed. Then, the confidence level of exclu- 
sion for planets with s and q is estimated as the fraction of bi- 
nary models not consistent with the best-fit single-lens model 
among all tested models with various values of a. For fit- 
ting binary models, it is required to produce many light curves 



with finite magnifications. We produce light curves by using 
the "map-making method" (iDong et alJl2006h . where a mag- 
nification map for a given s and q is constructed and light 
curves with various source trajectories are produced based on 
the map. In Figure [121 we present the obtained exclusion di- 
agrams for all analyzed events. Here we adopt a threshold of 
planet detection as Ax^^^ = Xs ~ Xp = 200, where Xp and 
represent the values for the best-fit planetary and single- 
lens models, respectively. For most events, the constraints 
on the excluded parameter space is not strong mainly due 
to the severe finite-source effect. However, the constraint is 
strong for MOA-2009-BLG-174 because of the small source 
size (p* ^ 0.002) and dense coverage of the peak. 

6. SUMMARY 

We provide integrated results of analysis for 14 high- 
magnification lensing events with lenses passing over the 
surface of source stars that have been detected since 
2004. Among them, 8 events are newly analyzed in this 
work. The newly analyzed events are MOA-2007-BLG- 
176, MOA-2007-BLG-233/OGLE-2007-BLG-302, MOA- 
2009-BLG- 174, MOA-20 1 0-BLG-436, MOA-20 1 1 -BLG- 
093, MOA-20 ll-BLG-274, OGLE-201 1-BLG-0990/MOA- 
201 l-BLG-300, and OGLE-201 1-BLG-l lOl/MOA-201 1- 
BLG-325. Information about the lenses and lensed stars ob- 
tained from the analysis is summarized as follows. 

1. For all newly analyzed events, we measure the linear 
limb-darkening coefficients of the surface brightness 
profile of the source stars. 

2. For all events with available CMDs of field stars, 
we measure the Einstein radii and the lens-source 
proper motions. Among them, 5 events (OGLE-2004- 
BLG-254, MOA-2007-BLG-176, MOA-2007-BLG- 
233/OGLE-2007-BLG-302, MOA-201 l-BLG-093, 
and MOA-201 l-BLG-274) are found to have Einstein 
radii less than 0.2 mas, making the lenses of the events 
candidates of very low-mass stars or brown dwarfs. 

3. The measured time scale Ie ^ 2.7 days combined with 
the small Einstein radius of ~ 0.08 mas of the event 
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Fig. 10. — Color-magnitude diagrams of neighboring stars in the fields of lensing events. In each panel, the circle represents the location of the lensed star and 
the triangle is the centroid of clump giants that is used as a reference for color and brightness calibration. For MOA-2009-BLG-174, the 'X' mark denotes the 
location of the blend. 



MOA-2011-BLG-274 suggests the possibility that the 
lens is a free-floating planet. 



4. For MOA-2009-BLG-174, we additionally measure the 
lens parallax and thus uniquely determine the physical 
parameters of the lens. The measured lens mass of ~' 
0.8 Mq is consistent with that of a star blended with the 
source, suggesting the possibility that the blend comes 
from the lens. 



We find no statistically significant planetary signals for 
any of the events analyzed in this work. However, it 
is still possible to place constraint on the range of the 
planetary separation and mass ratio. For this purpose, 
we provide exclusion diagrams showing the confidence 
levels of excluding the existence of a planet as a func- 
tion of the separation and mass ratio. 
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Fig. 12. — Exclusion diagrams of planets as a function of the star-planet separation (normalized in the Einstein radius) and the planet/star mass ratio. 



